MAGNETIC ACTIVITY IN STELLAR MERGER PRODUCTS 

Noam Soker^ and Romuald Tylenda^ 

ABSTRACT 

We study the expected X-ray luminosity of stellar merger products several 
years after merger. The X-ray emission is assumed to result from magnetic ac- 
tivity. The extended envelope of the merger product possesses a large convective 
region and it is expected to rotate fast. The rotation and convection might 
give rise to an efficient dynamo operation, therefore we expect strong magnetic 
activity. Using well known relations connecting magnetic activity and X-ray lu- 
minosity in other types of magnetically active stars, we estimate that the strong 
X-ray luminosity will start several years after merger, will reach a maximum of 
La; ~ 3 X 10^° erg s~^, and will slowly decline on a time scale of ~ 100 yr. We pre- 
dict that X-ray emission from V838 Mon which erupted in 2002 will be detected 
in 2008 with 20 hours of observation. 

Subject headings: stars: supergiants — stars: main sequence — stars: binary — 
stars: individual: V838 Mon — stars: magnetic activity — stars: merger 

1. INTRODUCTION 

The eruption of V838 Mon in 2002 (Brown 2002) and subsequent studies of its observed 
evolution (Munari et al. 2002; Kimeswenger et al. 2002; Crause et al. 2003; Kipper et al. 
2004; Tylenda 2005), as well as, of other similar objects, i.e. V4332 Sgr (Martini et al. 1999; 
Tylenda et al., 2005) and M31 RV (Mould et al. 1990) have led to suggestions that these 
observed events were likely to be due to stellar mergers (Soker & Tylenda 2003, Tylenda 
& Soker 2006). Soker & Tylenda (2006), who termed these events mergebursts, discuss the 
different channels to produce a mergeburst. 

For hours to months after merger the merger product is very luminous (e.g., Soker & 
Tylenda 2003; Bally & Zinnecker 2005; Tylenda 2005). For a grazing collision (namely, not 
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a head on collision) an extended envelope is inflated by the merging stars. Still on a longer 
time scale, the mass loss process, both mass loss rate and geometry, is strongly influenced 
by the merger event (e.g. Morris & Podsiadlowski 2006). On a much later time of hundreds 
of years and longer, after the merger products reaches equilibrium, the process can alter the 
evolution of the star on the HR diagram (e.g., Podsiadlowski et al. 1990), like the formation 
of blue stragglers (e.g. De Marco et al. 2005; Sills et al. 2005). 

In the present study we examine whether and when the merger product can become 
magnetically active, a process that might be observed in the X-ray and radio bands. 



2. THE EXPECTED MAGNETIC ACTIVITY 

2.1. The Rossby Number and X-ray Luminosity 

The parameter that best indicates the level of magnetic activity of main sequence stars 
(e.g., Pizzolato et al. 2003), pre- main sequence stars (e.g., Preibisch et al. 2005), and 
subgiants (or G giants; e.g., Gondoin 2005) is the Rossby number 

J-, -Prot Prot /-, \ 



where P^ot is the rotation period of the star, Tc = aHp/vc is the convection overturn time. Hp 
is the pressure scale hight, aHp is the mixing length, and is the velocity of the convective 
cells. In particular, the correlations of some properties of the magnetic activity in main 
sequence stars with the Rosbby number and the explanation of these in the frame of the 
auj dynamo model are well established (e.g., Brandenburg et al. 1988; Saar & Brandenburg 
1999). The subscript '6' in equation ([1]) indicates that the value of Tc is calculated at the 
bottom (inner boundary) of the envelope convective region, or just above it. In stars having 
a fully convective envelope, e.g., low mass main sequence stars, it is complicated to calculate 
Tc^f,. In that case one can deflne the global overturn time 

Tc-global = / , (2) 

jRb 

where i?* is the stellar radius and Rb is the radius at the bottom of the envelope convective 
region. Kim & Demarque (1996) flnd for main sequence stars the relation Tc-b — 0.5rc_giobai- 

We are interested in the x-ray emission resulting from magnetic activity. The magnetic 
flux on the surface of magnetically active main sequence stars is proportional to the X-ray 
luminosity (e.g., Pevtsov et al. 2003). The ratio of the X-ray luminosity to the bolometric 
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luminosity of main sequence stars has a general relation of 

^ = C,Ro-^; 0.15 <Ro< 10. (3) 

Lx/L]^o\ saturates at a values of ~ 10~^ for Ro < 0.15, while no activity is detected for 
Ro > 10 (Pizzolato et al. 2003). For main sequence stars ~ 10"^ (Pizzolato et al. 2003), 
while for G giants (subgiants) Cx ~ 10~^ (Gondoin 2005). YSOs are usually in the saturated 
regime, and show higher activity than main sequence stars with the same mass or bolometric 
luminosity (Preibisch et al. 2005). 



2.2. The Rossby Number in Inflated Merger Products 

Following Tylenda & Soker (2006) we assume that the merger remnant is composed of 
a more or less undisturbed pre-merger primary star of mass. Mi, and radius, Ri, surrounded 
be an envelope of mass, Menv, inflated up to an outer radius, i?env 

Merger products are expected to contract more or less along the Hayashi line (Tylenda 
2005; Tylenda & Soker 2006). However, they are different from young stellar objects (YSOs) 
contracting along the Hayashi line. An inflated merger product has a well defined and 
relaxed central region— the pre-merger primary star, while the contracting envelope contains 
a relatively small amount of mass. 

In that respect, the inflated merger remnants are more similar to late asymptotic giant 
branch (AGB) and post-AGB stars; both classes of objects share the following properties: 

1. Radius of tens to hundreds solar radii. 

2. Luminosity of ~ 3 x 10^ — 1O^L0. 

3. Cool envelope, T^s < 10^ K. 

4. Extended convection region in the envelope. To compensate for the low density in the 
expression for convective energy transport, the convective velocity must be large. 

5. A low mass envelope with a compact massive center: the stellar core in late AGB stars 
and post AGB stars, and the primary in inflated merger products. 

Based on these properties, we proceed as follows. To estimate the convective velocity Vc 
we use results of late AGB and post-AGB stars (Soker & Harpaz 1992; 1999). These results 
show that just below the photosphere, where the temperature is T ~ 10^ K, the convection 



-4- 



velocity is f c ~ 8 km s~^. In the stellar numerical code the convection velocity is limited by 
the isothermal sound speed, because for higher convection velocities the dissipation is large, 
and the convective cells rapidly slow down. The value of Vc stays at — 8 — 20 km s~^ in 
most of the envelope. We will therefore take Vc = 10 km s~^, and use the global convective 
overturn time as defined in equation ([2]). Using Tc-b — 0.5rc_giobai (Kim & Demarque 1996), 
and Rb ^ Rtfi we take for merger remnants 

r_ ^ 0.5^ ^ 40-^days. (4) 

A similar result is obtained if we consider, following Tylenda (2005), the envelope of the 
merger product to be an n = 3/2 polytrope, and calculate Tc at the middle of the envelope 
R = i?env/2. In an n = 3/2 envelope the pressure scale hight has its maximum value of 
Hp ~ i?env/10 at the middle of the envelope. Taking for the ratio of mixing length to pressure 
scale hight a = 1.86 (Kim & Demarque 1996) would give Tc-m — 15(i?env/lOO-R0) days. On 
the other hand, our estimate of Vc might be too large, with an underestimate of r^, as 
pre-main sequence stars have ~ 200 day (Preibisch et al. 2005). 

The inflated envelope of the merger remnant stores an angular momentum comparable 
to that of the pre-merger orbital motion of the secondary. For an = 3/2 polytropic envelope 
having i?env ^ Ri the moment of inertia can be approximated as / ^ O.llMenv-Renv We 
assume that after several dynamical time scales the convection in the envelope brings the 
envelope to a solid body rotation. Assuming that the secondary had a Keplerian velocity as 
it collided with the primary at radius -Ri and that the merger product envelope has a mass 
comparable to that of the secondary, we can estimate a rotation period of the envelope as 

Equivalently we can define a parameter t] being the ratio of the envelope rotation velocity 
to the Keplerian velocity fKcp(or Keplerian period Pkgp to rotation period) at Rem, namely 

,^(^) .0.9«V^ (6) 



yKepJ D \ R, 

The second equality uses equation ([5]). 

As it is clear from the above equations, when the remnant contracts, it spins-up. We 
assume that after it reaches a rotation velocity of some fraction rj^ax of its break-up (Ke- 
plerian) velocity mass loss keeps the value of r] unchanged. When it happens, the rotation 
period is 

— (^)"'fe)"©"^-- 
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The Rossby number (eq. [T]) for the inflated merger remnant can be obtained from 
equation (jl]) using equations ([5]) or ([7]), i.e. 



-1/2 /ox -1/2 

i?o(merger) ~ 3 



( Ren. \ (MA ^ (Rl 



if equation (Ej) gives rj < rj^a^ or 



1 / Ren. y M, 



1/2 / X -1/2 



otherwise. 



2.3. The X-Ray Luminosity of Inflated Merger Products 

As the merger products are somewhat similar to giant stars, we should take = 10~^ in 
equation ([3]) (Gondoin 2005). The operation of an au dynamo in the envelope of AGB stars 
that were spun-up by low mass companions spiraling inside their envelope was considered 
before (Nordhaus & Blackman 2006 and references therein). However, AGB stars that are 
expected to rotate very slowly and have large Rossby number Ro 3> 10 (Soker & Zoabi 2002), 
do amplify magnetic fields, as evidenced by polarization of maser emission in local regions 
around these stars (Szymczak 1998; Vlemmings 2005). In seems as if a dynamo based mainly 
on convection, and not on convection+rotation (the aVt dynamo model), can also amplify 
magnetic fields in giants (Soker & Zoabi 2002; Soker & Kastner 2003; Dorch 2004), but not 
as efficiently as the auj dynamo we appeal to here. Therefore, although our envelope model 
is similar to that of AGB stars, the dynamo model we use is much more efficient than that 
expected in AGB stars. By taking = 10~^ we might underestimate the X-ray luminosity 
of merger remnants. Using equation ([8]) or ([9]) in equation ([3]) with Cx = 10~^ we find the 
expected X-ray luminosity of the contracting envelope 

^ 4 X 10'" " erg (10) 

if equation ([HD gives r] < ri^s,x, or 

Lx ^ 1.2 X 10^° H-^Y ( "7^ V erg s^, (11) 



otherwise. 
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3. RESULTS FOR V838 Mon 

We can apply the general derivation of the previous section to predict the expected 
evolution of the X-ray luminosity of V838 Mon. 

As discussed in Tylenda (2005) the observed decline in flux of V838 Mon after its 
eruption can be well described by gravitational contraction of a low-mass inflated envelope 
sitting on top of an early B-type main sequence star. Assuming more recent determinations 
of the distance to V838 Mon giving a value of ~ 6 kpc (Sparks et al. 2007; Bond & Afsar 
2007) (compared to 8 kpc assumed in Tylenda 2005) the parameters of the model fitted to the 
observed decline become: Mi ~ TMq, Ri ~ 3.5i?0 and Mgnv — O.12M0. At the beginning of 
the contraction (August-September 2002) the envelope radius was i?env — 2OOOi?0. Using the 
same approach as in Tylenda (2005) and the above parameters, we can follow the contraction 
of the V838 Mon remnant to obtain the evolution of i?env and Lboi with time. This allows 
us to predict from the relations derived in Section [2] the evolution of the X-ray luminosity. 
The results are presented in Fig. [H The curves show the evolution of the X-ray luminosity 
with time and are labelled with the value of r/max used when caluculating the luminosity 
from equations (fTOj) and (ITT!) . The rising parts of the two upper curves correspond to the 
initial, constant angular momentum phase of the remnant contraction (77 < r^max)- The 
declining parts show the phase when the envelope is loosing angular momentum via mass 
loss so that the condition r] = rj^ax is kept. The dotted parts of the curves show the periods 
when the Rossby number is greater then 5. We expect that during this period the dynamo 
is less effective than assumed in our estimates so our results may overestimate the X-ray 
luminosity. 

With Chandra 23 hours of observations of the Orion Nebula, at a distance of 0.45 kpc, 
Feigelson et al. (2002) could detect sources with luminosity down to = 10^^ erg s~^. 
For a distance of 6 kpc to V838 Mon (Sparks et al. 2007; Bond & Afsar 2007), and with a 
similarly long observation, we expect to detect any emission if > 2 x lO^'' erg s~^. We 
conservatively took C^^ = 10~^ in equation ([3]), as appropriate for subgiants (Gondoin 2005) 
rather than = 10~^ as appropriate for main sequence stars (Pizzolato et al. 2003). More 
than that, bright pre main sequence stars with no accretion disk are X-ray brighter than 
those with disks (Preibisch et al. 2005). As V838 Mon does not have an accretion disk, it is 
quite possible that we underestimate the X-ray luminosity of merger products in equations 
flTU]) and flTT]) by up to an order of magnitude. Therefore, it is quite possible that 10 hours of 
XMM-Newton or Chandra observation could detect X-rays from V838 Mon at present and 
in coming years. 

V838 Mon was observed with Chandra for 6800 s a year after its outburst by Orio et al. 
(2003) who were able to put only an upper limit of Fx < 6.5 x 10~^'^erg cm~^ s~^. With 
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Fig. 1. — The expected evolution of the X-ray luminosity from the magnetically active 
V838 Mon remnant according to the merger model. The curves are labelled with the value 
of 77max; the maximum ratio of rotation velocity to Keplcrian velocity on the equator. Dotted 
parts indicate the time period when the Rossby number is Ro > 5, where the acu dynamo 
model is less efficient. 
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a distance of ~ 6 kpc this corresponds to Lx < 2.8 x lO^^ergs ^ which is well above our 
predictions. 

4. SUMMARY 

According to the stellar merger model of the V838 Mon outburst and similar merger 
products (which we term mergebursts), a large envelope is formed around the more massive 
of the two merging stars. The envelope then contracts on a thermal time scale. The merger 
remnant should become a fast rotator as it contracts. As the remnant contracts more or 
less along the Hayshi line, its envelope possesses a large convective region. The fast rotation 
and the envelope convection are the two ingredients required in the au dynamo model— a 
successful model for magnetic activity of main sequence stars, pre-main sequence stars, and 
subgiants. 

We applied the au model to contracting merger products by using the Rossby number 
(eq. [1]), and the relation between the Rossby number and X-ray luminosity known for 
magnetically active stars, scaled according to the expression for subgiant (or G giant) stars 
(eq. [3]). We also assumed that after the contracting product reaches some fraction r/max of its 
break-up (Keplerian) velocity, this ratio does not increase any more, because a stellar wind 
removes angular momentum from the envelope. Our final (and conservative) prediction for 
the X-ray luminosity of magnetically active merger products are given by equation (fTOl) for 
merger products before they reach our assumed maximum rotation rate, and by equation 
(fTTl) for merger products rotating at r^max- 

In section El we apply the results to our model of V838 Mon. The results are presented 
in Fig. [1] for three values of the assumed maximum rotating rate r^max, as marked near the 
lines. For too large Rossby numbers Ro > 10 (Pizzolato et al. 2003; we here take a stronger 
constraint of Ro ^ 5) of the au dynamo is not efficient any more. The dotted lines are 
the evolutionary stages where the expected Rossby number of V838 Mon is Ro > 5, and we 
expect no strong magnetic activity. 

[FTom Fig. [1] we learn the following. 

1. There is no magnetic activity at the first several years, and hence no X-ray emission is 
expected. The observation by Orio et al. (2003) was made a year after the outburst, 
when no magnetic activity and no X-ray emission is expected. 

2. For a reasonable values of maximum rotation rate 0.4 < r^max ^ 0.8, V838 Mon will 
reach a maximum activity at 6-8 years after outburst. The expected X-ray luminosity 
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then slowly declines. 

3. The X-ray luminosity in the coming years will be ~ 3 x 10^° erg s^^. At the distance 
of V838 Mon the expected X-ray flux is Fx ~ 6 x 10~^^ergcm~^ s~^. We estimate that 
with 100,000 seconds of observation this emission can be detected. 

We therefore highly encourage 100,000 seconds of X-ray observation of V838 Mon in 
2008. Even if no X-ray is detected, the results is of some importance, as it can strongly 
constrain models for V838 Mon, e.g., rules out accreting white dwarf. Orio et al. (2003) 
noted that their null detection rules our a symbiotic-like event to the V838 Mon outburst. 

We point out that the null detection of X-ray emission from two AGB stars (Kastner 
& Soker 2004) is not directly relevant to the case of V838 Mon. First, and most important, 
our prediction is based on the auo dynamo model, namely, the amplification of the magnetic 
field by the operation of both rotation and convection, which is known to be very efficient. 
On the other hand, predictions for AGB stars are based on the amplification of the magnetic 
field by convection alone (Soker & Zoabi 2002), which is thought to be much less efficient. 
Second, V838 Mon is an order of magnitude more massive than an upper AGB star. We 
predict the magnetic activity to take place when the radius, luminosity and temperature of 
V838 Mon is similar to that of upper AGB star. Due to the higher mass we expect the mass 
loss rate to be smaller, and the wind speed to be faster. Therefore, the column density to 
the expected X-ray emitting region will be much lower. 

Finally, the magnetic fields might be detected also in masers spots. Deguchi (2005) and 
Claussen (2005) report the detection of SiO maser around v838 Mon. We predict that if 
maser emission, in SiO, H2O, or OH, will be observed from 2007, some regions might show 
polarization indicating the presence of magnetic fields, similar to the case around AGB stars, 
e.g., Vlemmings et al. (2005). 

We thank an anonymous referee for useful comments. This research was supported in 
part by the Asher Fund for Space Research at the Technion, as well as, from the Polish State 
Committee for Scientific Research grant no. 2 P03D 002 25. 
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